GLOBAL DYNAMICS ABOVE THE GROUND STATE ENERGY 
FOR THE FOCUSING NONLINEAR KLEIN-GORDON EQUATION 

K. NAKANISHI AND W. SCHLAG 

Abstract. The analysis of global dynamics of nonlinear dispersive equations 
has a long history starting from small solutions. In this paper we study the 
r»o ■ focusing, cubic, nonlinear Klein-Gordon equation in M 3 with large radial data in 

the energy space. This equation admits a unique positive stationary solution Q, 
called the ground state. In 1975 Payne and Sattinger showed that solutions u{t) 
with energy E[u,u] strictly below that of the ground state are divided into two 
classes, depending on a suitable functional K(u): If K (u) < 0, then one has finite 
time blow-up, if K(u) > global existence; moreover, these sets are invariant 
Q^ ■ under the flow. Recently, Ibrahim, Masmoudi and the first author [22] improved 

-^ this result by establishing scattering to zero for K[u] > by means of a variant 

of the Kenig-Merle method [25], [26]. In this paper we go slightly beyond the 
ground state energy and we give a complete description of the evolution in that 
case. For example, in a small neighborhood of Q one encounters the following 
trichotomy: on one side of a center-stable manifold one has finite-time blow-up 
for t > 0, on the other side scattering to zero, and on the manifold itself one has 
scattering to Q, both as t — > +oo. In total, the class of data with energy at most 
slightly above that of Q is divided into nine disjoint nonempty sets each displaying 
different asymptotic behavior as t-s- ±oo, which includes solutions blowing up in 
one time direction and scattering to zero on the other. The analogue of the 
solutions found by Duyckaerts, Merle [13], [14] for the energy critical wave and 
Schrodinger equations appear here as the unique one-dimensional stable/unstable 
manifolds approaching ±Q exponentially as t — > oo or t — > — ex), respectively. The 
main technical ingredient in our proof is a "one-pass" theorem which excludes the 
existence of (almost) homoclinic orbits between Q (as well as — Q) and (almost) 
heteroclinic orbits connecting Q with —Q. In a companion paper [31] we establish 
analogous properties for the NLS equation. 
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1. Introduction 

In this paper we study the global behavior of general solutions to the nonlinear 
Klein-Gordon equation (NLKG) with the focusing cubic nonlinearity on IR 3 , i.e., 

u-Au + u = u 3 , u(t,x) : M 1+3 -fR, (1.1) 

which conserves the energy 

f U'/l 2 -|_ IV7/I 2 -I- I7/I 2 I7/I 4 
E (u) := / \ M +|V o M| +N - qL] dx. (1.2) 



We regard /J 1 (1R 3 ) x L 2 (R 3 ) as the phase space for this infinite dimensional Hamil- 
tonian system. In other words, we write the solutions as 

u(t) := (u(t),u(t)) en:= Hl d x L r 2 ad - (1.3) 

There exists a vast literature on the wellposedness theory for this equation in the 
energy space since Jorgens [23] and Segal [35], as well as the scattering theory 
for small data for the focusing nonlinearity as in (1.1) and for large data for the 
defocusing equation; see Brenner [7], [8], Ginibre, Velo [16], [17], Morawetz, Strauss 
[30], and Pecher [33]. In this paper, the scattering of a solution u to a static state 
tp refers to the following asymptotic behavior: there exists a solution v of the free 
Klein-Gordon equation such that 

\\u(t)-0-v(t)\\ H ^O, (1.4) 

as either t — >■ 00 or t — > —00 (depending on the context). See also Strauss [37] 
and [22] for a review of Strichartz estimates and wellposedness, as well as scattering 
in this setting. 

It is well-known [36], [5], [11] that there exists a unique radial positive ground 
state Q(x), solving the static equation 

-AQ + Q = Q 3 , (1.5) 

with the least energy 

E(Q) = m := / [ w±w _ jgr dx > o> (16) 

among the static solutions, and that the solutions u below the ground energy 

E(u) < E{Q) (1.7) 
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are split into two classes by the functional 

K (u) := f [\Vu\ 2 + \u\ 2 - M 4 ] dx = d x \ x=1 J(\u), (11 



(1) If Ko(u(0)) > then the solution u exists globally on t G R. 

(2) If K (u(0)) < then the solution u blows up both in ±t > 0. 

In this paper, blow-up of a solution u means that it cannot be extended beyond a 
finite time T* in the energy space %. Then the wellposedness theory implies that 
||u(£)||tt — > oo as t approaches T*. 

The above result of dichotomy is essentially due to Payne, Sattinger [32], which 
of course implies that the ground state Q is unstable. For a general theory of 
orbital stability vs. instability of solitary wave solutions for equations of this type, 
see Grillakis, Shatah, Strauss [18], [19]. 

Recently, Kenig, Merle [25, 26] considered the focusing nonlinear Schrodinger 
equation (NLS) and the nonlinear wave equations with the energy-critical nonlin- 
earity in three dimensions (as well as others) and obtained a stronger version of the 
dichotomy, by adding scattering to zero in the global existence scenario. It has been 
extended to several similar equations, including the cubic focusing NLS equation by 
Holmer, Roudenko [21], and NLKG such as (1.1) in [22]. 

On the other hand, the second author [34], followed by Beceanu [4], constructed 
center-stable manifolds with finite codimensions around the ground state in the 
case of the focusing cubic NLS equation in three dimensions. A construction of 
such manifolds for all L 2 supercritical NLS equations in one dimension was given by 
Krieger and the second author [27], and for the radial critical wave equation in M 3 
this was done in [28]. Stable, unstable and center manifolds have been known to 
play an important role in hyperbolic dynamics for a long time, see the classical work 
by Hirsch, Pugh, Shub [20] as well as many other works since then, such as Bates, 
Jones [2], [3], Vanderbauwhede [38], and Li, Wiggins [29]. While the more ODE 
oriented approach in Bates, Jones to PDEs such as (1.1) does construct center-stable 
and unstable local manifolds near an unstable equilibrium in the energy space, no 
statement can be made about the trajectories once they leave a small neighborhood 
of the equilibrium. In other words, they are local-in-time results which do not 
involve any dispersive analysis. In contrast, the emphasis in [34], [27], [28], [4] lies 
with the asymptotic behavior of the solutions as t — > oo and scattering to a suitable 
"soliton" is established in each of these works for data lying on a center-stable 
manifold modulo the group actions which preserve the equation. In [4] both points 
of view are united by establishing the existence of the manifold Ai from [34] in the 
energy class (in fact, in the critical H 2 topology) for the cubic NLS equation, with 
the added feature of global-in-time invariance of Ai as t — > 00, as well as scattering 
in the energy class to Q modulo the symmetries for solutions starting on Ai. 

Note that those solutions scattering to Q must have energy above the ground 
state, and so the solution sets of Payne, Sattinger on the one hand, and those on 
the center-stable manifolds on the other hand, are necessarily disjoint from each 
other. More recently, Duyckaerts, Merle [13, 14] investigated the solutions on the 
threshold energy E{u) = E(Q) for energy critical equations, and proved that there 
are exactly two new solutions modulo symmetry; scattering to Q as t — > —00, while 
either scattering to or blowing up as t — > 00. They can be regarded as minimal 
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solutions on the manifold Ai. This result is also extended to the cubic NLS by 
Duyckaerts, Roudenko [15]. In this paper we exhibit these solutions as the unique 
stable/unstable manifolds associated with ±Q, see [3] for the definition of these 
objects. 

However, all of these works describe only part of the dynamics for energies near 
that of the ground state. A natural question to ask is whether the solutions near the 
soliton are separated by M. into a region of scattering to zero, and the other region 
of blow-up. This is partially motivated by the study of "critical phenomena" in the 
physics literature, see for example Choptuik [9], Choptuik, Chmaj, Bizon [10], and 
Bizon, Chmaj, Tabor [6]. For the |w| 5 wave equation in R 3 , blow-up for all data on 
one side of the tangent space to the center-stable manifold constructed in [28] was 
shown by Karageorgis, Strauss [24]. 

Our goal in this paper is to give a complete picture of the dynamics of (1.1) for 
radial energy data with energy slightly larger than that of Q: 

H £ :={ueU\E(u)<E(Q)+e 2 }. (1.9) 

Note that the only symmetry in "H is u \— > —u in this setting. 

Theorem 1.1. Consider all solutions of NLKG (1.1) with initial data u(0) € T-t £ for 
some small e > 0. We prove that the solution set is decomposed into nine non-empty 
sets characterized as 

(1) Scattering to for both t — > ±oo ; 

(2) Finite time blow-up on both sides ±t > ; 

(3) Scattering to as t — > oo and finite time blow-up in t < 0, 

(4) Finite time blow-up in t > and scattering to as t — > — oo ; 

(5) Trapped by ±Q for t — ¥ oo and scattering to as t — >■ — oo ; 

(6) Scattering to as t —¥ oo and trapped by ±Q as t — > — oo ; 

(7) Trapped by ±Q for t — ¥ oo and finite time blow-up in t < 0, 

(8) Finite time blow-up in t > and trapped by ±Q as t — >■ — oo ; 

(9) Trapped by ±Q as t — >■ ±oo ; 

where "trapped by ±Q" means that the solution stays in a 0(e) neighborhood of ±Q 
forever after some time (or before some time). The initial data sets for (l)-(4), 
respectively, are open. 

The striking difference from the Kenig-Merle or Duyckaerts-Merle type results 
is the existence of solutions which blow up for t < and scatter for t — > +oo, 
or vice versa. It also implies that the initial data set for the forward scattering 
(1) U (3) U (6) (or backward scattering) is unbounded in H; in fact, it contains a 
curve connecting zero to infinity in %. The number "nine" simply means that all 
possible combinations of scattering to zero/scattering to ±Q/finite time blow-up are 
allowed as t — > ±oo. Each of these are in fact realized by infinitely many solutions. 

The simple dichotomy in terms of Kq is no longer available, so it is not easy to 
predict the global dynamics for a given initial data, as in the case of Payne-Sattinger 
or Kenig-Merle below the ground state. However, we can give some description, in 
a more dynamical way, combining the hyperbolic structure around the ground state 
with the global variational structure. Let L + denote the linearized operator around 
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Q, and let p be its ground state: 

L+ = -A + 1 - 3Q 2 , L + p=-k 2 p, k>0, p(x)>0, \\p\\ L 2 = l. (1.10) 

It is well known that L + has only one negative eigenvalue. We can define a nonlinear 
distance function (Iq(u) : "H e — > [0, oo) continuous such that 



d Q (u) ~inf ||u=fQ||«, 



;i.n) 



d Q («)<i => rf Q (M) = E(M)-E(g) + A; 2 |( MT g|p)| 2 ; 

where ±Q is chosen to be the closest to w. Let 

B(±Q) := {ueH £ \ d Q {u) < 2[E(u) - E(Q)}} (1.12) 

be a pair of small balls around ±Q. Outside of these balls, we can define a continuous 
sign function & : "H £ \ B(±Q) — > {±1} such that 

d Q (u) < S x =► 6(u) = -sign(w ^ Q|p), 
dQ(u) > 5 5 =^> 6(w) = sign^oH, 

for some 8x > Ss > 2s, with the convention that signO = +1. We emphasize that 
& is different 1 from signfQ) in the region close to ±Q. 

Theorem 1.2. There are small e > and R > 2e with the following property. If u 
is a solution of NLKG (1.1) in H 6 , defined on an interval I, satisfying 

d Q (tf(n)) < d Q (u(r 2 )) < R, u(t 2 ) £ B(±Q), (1.14) 

for some T\ < r 2 G /, then we have dq(u(t)) > ^(w^)) for all t > r 2 in I . In 
particular, u remains outside of B(±Q) with a fixed sign &(u(t)) G {±1} fort > r 2 . 
If &(u(t)) — +1, then u scatters to as t — > oo. If &{u(t)) = —1, then u blows 
up in finite time after r 2 . Conversely, if a solution u in l-i e scatters as t — > oo ; 
then dq{u{t)) > R, u(t) $_ B(±Q) and &(u(t)) = +1 for large t. If it blows up as 
t ->■ T - 0, then d Q (u(t)) > R, u(t) £ B(±Q) and &(u(t)) = -1 for t < T close 
to T. If it is trapped by ±Q as t —¥ oo ; then it stays in B(±Q) for large t. 

In other words, every solution can enter and exit B(±Q) at most once, and the 
sign of &(u) is constant while it is away from B(±Q). Hence, if a solution does 
enter B(±Q) and exits, then its fate is determined by the sign of & at the entrance 
and the exit, respectively, for t — >■ — oo and t — > +oo, which are not necessarily the 
same (cf. cases (3) and (4) of Theorem 1.1). It also implies that the set (9) in the 
previous theorem is bounded in the energy space %. 

In fact, we can give a more precise description of the exiting dynamics. Indeed, 
in the setting (1.14), dq(u(t)) is monotonically and exponentially growing for t > r 2 
until it reaches a larger number Sx (see Lemma 4.2). 

One can also give a more detailed description of the dynamics of those solutions 
which are trapped by ±Q, by means of the following spectral gap property of the 



& might be the same as sign_Ko in H. 6 \ B(±Q), but it depends on the coefficient 2 in the 
definition of B(±Q). In fact we can make it closer to 1 as cZq(m) — > 0, and in that region we can 
show that — sign(u =p Q\p) and signifo(u) are indeed opposite at some points. 
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linearized operator L + = —A + 1 — 3Q 2 : 

L + has no eigenvalues in (0, 1] and no resonance 2 at the threshold 1 (1.15) 

It was verified by means of the Birman-Schwinger theorem and a numerical com- 
putation of the five largest eigenvalues of a suitably discretized Birman-Schwinger 
operator in [12] by Demanet and the second author. Using this gap property one 
obtains the following refinement of Theorem 1.1. 

Theorem 1.3. In the statement of Theorem 1.1 one can replace "trapped by ±Q" 
with "scattering to ±Q ". The sets (5) U (7) U (9) and (6) U (8) U (9) are smooth 
codimension one manifolds in the (radial) phase space T-L, and they are the 3 center- 
stable manifold, resp. the center-unstable manifold, around ±Q. Similarly, (9) is a 
smooth manifold of codimension 2, namely the center manifold. 

It seems natural to the authors to separate the above theorem from the previ- 
ous ones, as it is know from [12] that (1.15) fails if one lowers the power 3 on the 
nonlinearity slightly, say to power < 2.8. On the other hand, our argument for 
Theorems 1.1 and 1.2 is quite general, in particular it does not use (1.15), or more 
precisely, any dispersive property of the linearized operator. In fact, it is straight- 
forward to extend them to all L 2 super-critical and H 1 subcritical powers and all 
space dimensions, i.e. 

u-Au + u = u p , l + A/d<p< l + 4/(d-2), u : R 1+d ->■ R. (1.16) 

Finally, from the above theorems, we can easily deduce a Duyckaerts-Merle-type 
result on the energy threshold. 

Theorem 1.4. Consider the limiting case e — > in Theorem 1.1, i.e., all the radial 
solutions satisfying E{u) < E(Q). Then the sets (3) and (4) vanish, while the sets 
(5) — (9) are characterized, with some special solutions W±, as follows 

(5) = {±W-(t - t ) | t G R}, (6) = {±W.(-t - t ) | t G R}, 

(7) = {±W + {t - t ) | t G R}, (8) = {±W + {-t - t ) | t G R}, (1.17) 

(9) = {±Q(t - t ) | t G R}. 

The solutions W±(t) converge exponentially to Q as t — >■ oo. 

As mentioned before (5) U (7) U (9) is the stable manifold, and (6) U (8) U (9) the 
unstable manifold, associated with ±Q, cf. [3]. The rest of the paper is organized as 
follows. In Section 2, we recall several known facts about the ground state and the 
linearized operator around it. In Section 3 we construct the center-stable manifold 
around the ground states, under the aforementioned spectral condition (1.15). Next 
to Theorem 1.3 this is the only part where we require the gap property. Section 4 is 
the most important and novel part of this paper, where we prove a part of Theorem 
1.2 that every solution can enter and exit a small neighborhood of the ground states 



In our radial case, it is easy to preclude by means of analytical arguments eigenvalues at the 
threshold 1, but we include it here, because such an eigenvalue might exist in the nonradial case, 
which is indeed treated by [12]. 

Since center manifolds are in general not unique it might be more precise to say "a center-stable 
manifold" here, but we ignore this issue. In fact, our manifolds arc naturally unique for the global 
characterization in Theorem 1.1. 



(2.3) 
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±Q at most one time. We call it "one-pass theorem". One then immediately obtains 
the blow-up part of the dynamical classification by the classical Payne, Sattinger 
argument. In Section 5, we then prove the scattering part, by complementing the 
Kenig, Merle argument with the one-pass theorem. Finally in Section 6, we describe 
the global dynamics and its classification as in Theorem 1.1 and Theorem 1.3, using 
some simple topological arguments and the one-pass theorem. In the appendix we 
give a table of notation for frequently used symbols. 

2. The ground state 

Here we recall several known properties of the ground state Q. First we consider 
its variational character with respect to the scaling symmetry. For our purpose, it 
suffices to consider the L° and the L 2 invariant scalings: 

<p(x) H- (Po(x) := vip(x), ip(x) H- <f 2 (x) := u 3/2 ip(ux). (2.1) 

Let D and D 2 , respectively, be the generators of these symmetries, viz. 

x^7 ~\~ ^7 x 
D ip(x) := ip(x) = d u \ u=1 ipQ, D 2 (p(x) := (p(x) = d u \ u=1 ip 2 , (2.2) 

and let Kq and K 2 be the derivatives of J with respect to these scalings: 
K (<p):= [ [\V ( p\ 2 + \<p\ 2 -\ ( p\*]dx = d v \ v=1 J(<p v ) = (J'( ( p)\D 0( p), 

</R 3 

K 2 ( V ): = [ [\V<p\ 2 -j\<p\']dx = d v \^ 1 J(<fi) = (J'(<p)\D 2 <p), 

where J'(<p) denotes the Frechet derivative 

J'(ip) = -Ay? + <p - if 3 . (2.4) 

Define positive functionals G s for s = 0, 2 by 

GM --=J-^ = \y\\h, G 2 (<p) ■■=J-^ = \\\V<p\\h + \\\<p\\h- ( 2 - 5 ) 
Lemma 2.1. For s = 0,2 we have 

J(Q) = inf{ J(<p) \0^<peH\ K S ( V ) = 0} 
= inf{G^) \0^<peH\ K S ( V ) < 0}, 
and these infima are achieved uniquely by the ground states ±Q. 

Proof. If K s (ip) < then K s (\*tp) = for some A* G (0,1), whereas G s (\*ip) < 
G s (ip). Since moreover K s (<p) = implies J((p) = G s (ip), the two infima are equal. 
Since J'(Q) = implies K S (Q) = 0, the infima are no larger than J(Q). 

To obtain a minimizer, let {<p n }n>i C H]^ \ {0} be a minimizing sequence such 
that (the Schwartz symmetrization allows us to restrict them to the radial functions) 

K s (ip n ) = 0, J(cp n ) -»■ m, (2.7) 

where m denotes the right hand side of (2.6). Since G s ((p n ) = J(<p n ) is bounded, (p n 
is bounded in H 1 . After extraction of a subsequence, it converges weakly to some 
tpoo in H 1 , and in the strong sense in L 4 , by the radial symmetry. Thus X s (</?oo) ^ 0, 
J^oo) < J{Q) an d Gg^ipoo) < G S (Q). If y?oo = then the strong convergence in 



(2.6) 
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L 4 together with K s ((p n ) — > implies that tp n — > strongly in H 1 . Notice that 
K s (ip) > for < Hv^lli/ 1 <C 1, due to the interpolation inequality 

Iklll* < Ikll^llV^lia. (2.8) 

Hence K s (ip n ) > for large n, a contradiction. Thus we obtain a nonzero minimizer 
(Poo, and so K^ip^) = 0, which implies that ip n — > ipoo strongly in H 1 . 

The constrained minimization implies that for some Lagrange multiplier \x G IR, 

j'^oo) = IlKfaoo), K s {Voo) = 0, (2.9) 

Multiplying the Euler-Lagrange equation with D s (p and integrating by parts yields 

= K s { Voo ) = rtK'.MD.^)) = » x {H 0{ ^\ " f^f I" = °|' (2.10) 

\2K 2 {(p 0O ) - f|boo|ll4 (s = 2), 

which implies /i = 0, so ip^ is a ground state. □ 

Next, we record the following observation from [22]. 

Lemma 2.2. For any (p G H 1 such that J(tp) < J{Q), one has either tp = 0, or 
K s (<p) > for both s = 0,2, or K s (<p) < for both s = 0, 2. 

Proof. In the proof of the above lemma, we already know that is surrounded by 
the open set K s (ip) > for both s, and K s (ip>) = is prohibited except for by 
J {(f) < J{Q), due to the above minimization property of Q. Then it is enough to 
show that each set /C+ := {K s > 0, J < J(Q)} is contractible to {0}, which implies 
that it is connected and so cannot be divided by the disjoint open sets /C^l s and 
/Cr_ s = {K 2 ~ s < 0, J < J(Q)}. The contraction is given by ip u s with v : 1 — > 0. By 
definition d v J{ip u s ) = K s (tp v s )/v, and so J(<p^) decreases together with u, as long as 
K s {ip> v s ) > 0, which is preserved as long as J(<p") < J(Q)- Hence, <p> v s stays in /C+ for 
1 > v > 0. When s = 2, it converges to as v — \ +0 only in H 1 fl L 4 , but since /C 2 
are open in that topology this is sufficient: once (p 2 gets in a small ball around in 
K-2, one can change to the other scaling (p i— >■ i/</? to send it to in if 1 . □ 

Next we recall the spectral properties of Q and the linearized operator L + defined 
in (1.10). Decomposing the solution of (1.1) in the form 

u = Q + v, (2.11) 

we obtain the equation of the remainder 

v + L + v = N(v), N(v) := (v + Q) 3 - Q 3 - 3Q 2 v = 3Qv 2 + v 3 . (2.12) 

The energy functionals are expanded correspondingly 

J(Q + v) = J(Q) + (J'(Q)\v) + (J"(Q)v\v}/2 + 0(\\v\\ 3 Hl ) 

= J(Q) + (L + v\v)/2 + 0(\\v\\ 3 Hl ), (2.13) 

K S (Q + v) = (K' s (Q)\v) + (K':(Q)v\v}/2 + 0(\\v\\ 3 H1 ). 
In particular 

MQ + v) = 2{-AQ + Q- 2Q 3 \v) + ((-A + 1 - 6Q 2 )v\v) + 0(\\v\\ 3 Hl ] 
= -2(Q 3 \v) + ((L+ - 3Q 2 )v\v) + Odlvlllfi). 



(2.14) 
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Lemma 2.3. As an operator in L 2 ad , L + has only one negative eigenvalue, which 
is non- degenerate, and no eigenvalue at or in the continuous spectrum [l,oo). 

Proof. L + has at least one negative eigenvalue because 

(L + Q\Q) = -Z\\Q\\\ < 0, (2.15) 

and at most one, because 

(Q 3 \v) = => (L + v\v)>0. (2.16) 

To see this, suppose that / G H 1 satisfies (Q 3 \f) = and (L + f\f) = — 1, and let 
v = £ Q + 5f for small e, 5 G R. Then 

MQ + v) = -2e\\Q\\ 4 L4 - 5 2 (l + (3Q 2 f\f)) + 0[e 2 + s5 + 5 3 ), (2.17) 

so there exists e = 0(S 2 ) such that Kq(Q + v) — 0. On the other hand 

J(Q + v) = J(Q) -5 2 + 0(5 3 ) < J(Q), (2.18) 

which contradicts the minimizing property Lemma 2.1. 

Next, if 7^ / G L^ad solves L+f = 0, then / _L Q 3 , Q, because L + Q = -2Q 3 
and L + (rd r + 1)Q = — 2Q- By the Sturm-Liouville theory, / can change the sign 
only once, say at r = r > 0, so does Q 3 — Q(r ) 2 Q, since Q(r) is decreasing. Then 
(Q 3 — Q(ro) 2 Q)f is non-zero with a definite sign, contradicting / _L Q 3 , Q. 

The absence of embedded eigenvalue is standard, and follows also from the as- 
ymptotic equation for any eigenfunction, viz. 

L + f = \ 2 f =► (d 2 + \ 2 )(rf) = -3Q 2 (rf) < e~ 2r , (2.19) 

by the exponential decay of Q. □ 

Let 

P + :=l-p(p\ (2.20) 

be the orthogonal projection for p. Then the above lemma implies (L + v\P + v) > 
\\v \\ 2 L 2, and so for any v G P + (H l ) and ^ G (0, 1] sufficiently small, 

(L + v\v) > (1 - 0)||v||2 + ^llvll^ - 3£(Q 2 ^> > || V ||^, (2.21) 

hence (L + v\v) ~ ||^||^-i on P + (H 1 ). 

The above property of L + is sufficient for the analysis of dynamics away from Q, 
but for that of solutions staying forever around the ground state, we require the 
following property of L + : 

L + has no eigenvalue in (0, 1] and no resonance at the threshold 1. (2.22) 

Both parts of this statement have been verified by Demanet and the second author 
in [12] via numerics. Their approach is based on the Birman-Schwinger theorem 
which equates the number of eigenvalues and resonances < 1 (the threshold) of 
L + , counted with multiplicity, to the number of eigenvalues > 1 of the self-adjoint, 
positive, compact operator K + := 3Q(— A)~ X Q, again counted with multiplicity. In 
fact, [12] by means of a numerical computation finds that K + has precisely four 
eigenvalues greater than 1 (corresponding to the negative ground state plus the zero 
eigenvalue of multiplicity three: L + VQ = 0), whereas the fifth largest eigenvalue 
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was calculated to be < 0.98 with an estimated 8 to 9 digits of accuracy behind the 
comma (to be precise, A5 = 0.97039244...). By the Birman-Schwinger theorem, 
this verifies (2.22). 

3. Center-stable manifold 

In this section, we investigate the solutions staying around the unstable ground 
state, constructing the center-stable manifold, while the center-unstable manifold is 
obtained by reversing the time. The situation is much simpler than for NLS [34], 
because the only symmetry present (reflection and time translation) fixes the ground 
state and the linearized operator is scalar self-adjoint. We decompose any solution 
u simply by putting 

u{t) = Q + v(t), v(t) = \(t)p + -y, 7^P, (3.1) 

which is obviously unique. We then obtain the equations of (A, 7) 6lx P + {H 1 ) 

'\-k 2 \ = P p N{v) =:N p (v), 



j + u 2 -f = P + N(v)=:N c (v), u:=y/P+L^. 



(3.2) 



We look for a forward global solution which grows at most polynomially, to which 
it is equivalent to remove the growing mode e kt . From the integral equation of A 
one extracts the growing mode 



1 



A(t) = cosh(A;t)A(0) + - sinh(fc£)A(0) + - / sinh(A;(t - s))N p (v)(s) ds 

t kJo (3.3) 

A(0) + -i-A(O) + \ f e- ks N p (v)(s) ds +■■■ 
k k Jo 



e kt 



where the omitted terms are exponentially decaying. Hence, the necessary and 
sufficient stability condition is 

/•oo 

A(0) = -jfeA(O) - / e- ks N p (v)(s)ds. (3.4) 

Jo 

Under this condition, the integral equation for X(t) is reduced to 

r 1 /-oo 

\{t) = e~ kt A(0) + — J e- ks N p (v)(s)ds 



l^ ,' e- k ^N p (v)(s)ds, (3.5) 







while the integral equation for 7 is 

7(t) = cos(wt)7(0) + - sin(wt)7(0) + - / sin(a;(t - s))N c (v)(s) ds. (3.6) 



1 „, 1 " 



The coupled equations (3.5)-(3.6) can be solved by iteration, using the Strichartz 
estimate for e ltuJ , for any small initial data (A(0), 7(0)) G R x P + ("H). The linearized 
energy norm of v around Q is denoted by 



Me := y/[k*(v\p) 2 + \\uP+v\\h + ll*IIL.]/2 , N 

/ : (3 ' 7) 

= V / [P|A|2 + |A|2]/2+||7|||. 
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By Lemma 2.3 we have 

\W\\ 2 E^\\v\\ 2 n =\\v\\ 2 H1 + \\v\\ 2 L2 . (3.8) 

Recall that % is the radial energy space. 

Proposition 3.1. Assume that (2.22) holds. Then there are v > and C > 1 with 
the following property: For any given A(0) G K, 7(0) G P + (7i) satisfying 

E := k 2 \\(0)\ 2 + \\«n(0)\\h + Il7(0)||| 2 < *A (3.9) 

t/iere exists a unique solution u of NLKG (1.1) on < £ < 00 satisfying 

n(0) = Q + A(0)p + 7 (0), P+n(0) = 7 (0), (3.10) 

|(p|n(£) + fcu(t))| ^ ^0 /or a// £ > 0, and 

\\u{t) - Q\\ 2 H < CE (0<V£<oo). (3.11) 

T/ie dependence of u on (A(0),7(0)) is smooth in L°°(0, 00; %). In addition, there 
exists a unique free Klein-Gordon solution 7^ such that 

\\{t)\ + |A(t)| + ||7(t) - 7oo(*)||« -> 0, (3.12) 

as t — > 00. In particular, we have E{u) = J(Q) + 1 1 Too 1 1 ^ /2 . 

Conversely, any solution u of NLKG satisfying (3.11) with E < u 2 /C must be 
given in this way, which is uniquely determined by (A(0),7(0)) ; and by (A(0),7oo(0)). 

An alternate, but essentially equivalent, formulation of this result reads as follows. 

Proposition 3.2. Assume that (2.22) holds. Then there exists v > small and a 
smooth graph Ai in B U (Q) cH so that Ai is tangent to 

T Q M = {(u , Ul ) en\ (ku + Ul \p) =0} (3.13) 

in the sense that 

sup dist(x, T Q M) <5 2 V <5 <v 

x&dBs(Q) 

and so that any data (uq,ui) G M. lead to global evolutions of (1.1) of the form 
u = Q + v where v scatters to a free Klein- Gordon solution in H. Moreover, no 
solution can stay off Ai and inside B V (Q) for all t > 0, and, At is invariant under 
the flow for all t > 0. 

Proof of Proposition 3.1. For the existence, we solve (3.5)-(3.6) by iteration using 
the norm 

||(A, 7 )|U := ||A||LinL~ ( o,oo) + H7lU(o,oc), St := L 2 L 6 X n L?Hl (3.14) 

The Strichartz estimate for the free Klein-Gordon equation gives us 

\\u\\st(o,T) < \\u(0)\\n + \\(d 2 - A + l)u\\ LlL% ^ T) . (3.15) 

Under the hypothesis (2.22), the operator L + satisfies the conditions of Yajima's 
W k,p boundedness theorem for the wave operators [39], so that we can conclude 
that (3.15) applies to 7 as well: 

ll7||st(o,T) < ll7(0)|k+||(«9 t 2 - W 2 )7||L^(o,r), (3.16) 
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provided that 7(0) and (d 2 — oo 2 ) r y are orthogonal to p. 

For the solution A of (3.5), with A(0) uniquely determined by (3.4), we can esti- 
mate the norm by simple integration in t: 

||A|UinL~(o,oo) < fc _1 |A(0)| + fe-lJV^IUi^oo) < |A(0)| + \\N(v)\\ L}Lmoo) , (3.17) 

and for the solution 7 of (3.6) by the above Strichartz estimate 

||7lko,oo) < Il7(0)|k + ||^(u)||li^(o,oo)- (3.18) 

The nonlinearity N(v) is bounded in L\L 2 X by 

\\Qv 2 \\l\li + \\A\l\l% < \H% L% {\\v\\ LTL% + \\Q\\l T l%), (3.19) 

where the norm of v is bounded by 

\Ml%l% < ||A|| L i n L t ~ + I|7||l?l« (1 < P < 00). (3.20) 

Gathering them, we obtain 

ll(A,7)IU < |A(0)| + ||7(0)|k + ll(A,7)llx + ll(A,7)llx- (3.21) 

Applying these estimates to the iteration sequence, we obtain a unique fixed point 
of (3.5)-(3.6) for any given small (A(0),7(0)). It is straightforward to see that 
u :— Q + X(t)p + 7(t) solves NLKG on < t < 00, satisfying 

ll^-Qlk<ll(A,7)Hx<|A(0)| + ||7(0)|k, (3.22) 

with smooth dependence on the data. The bound on (p\u+ ku) = X + kX follows 
by using the equation once again. 

Moreover, the asymptotic profile of 7 is given by 

1 1 f°° 

700(f) = cos(wt)7(0) + - sm(ut)j(0) + - / sin(w(* - s))N c (v)(s) ds, (3.23) 

w ■ u Jo 

with the convergence property 

||A|| l -(t,oo) + ||7-7oo|U t ~w(T,oo) ^ II^(^)I|lil2 (T)00 ) ^0 (T -»■ oo). (3.24) 
The iteration can be solved with a given 7^ and the equation of 7 now reads 

1 /•* 

l(t) = 7oo (t) + - / sin(a;(t - s))N c (v)(s) ds, (3.25) 



where the estimates are essentially the same. 7oo(^) can be further replaced with a 
free Klein-Gordon solution by the linear scattering for L + . 

The uniqueness part requires some more work, since a priori we do not know if the 
solution is in the space X, globally in time. Let u be a solution on [0, 00) satisfying 
(3.11). Since it is bounded in the energy space, we can easily see that N p (t) is 
bounded. Therefore, it has to satisfy (3.4), and the reduced integral equation (3.5) 
as well as (3.6) for all < t < 00. To see that A G L](0, 00) and 7 G L 2 L®(0, 00), 
consider the norm 

||(A,7)|U T := ||A|| L i(o,T) + ||A||l°°(t,oo) + ||7l|st(o,T), (3.26) 
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for T > 0. The energy bound implies that ||(A, 7)||x T < oo for all T > 0, but we 
require a uniform bound. From the integral equations one concludes that 



(3.27) 



||'M|,L 1 (0,T)nL° o (T,oo) ^S \\Np\\L 1 (0,T)nL°°(T,oo) 

< V + ^(||A|| L i(o,T)nL°°(T,oo) + ll7lUt(0,T)), 

and using the Strichartz estimate, one further has 

\\l\\st(0,T) < V + ^(||A||z,l(o,T) + \\l\\st(0,T)), (3.28) 

Thus we obtain by Fatou 

||(A,7)|| x <liminf||(A, 7 )|U T <v, (3.29) 

and the contraction mapping principle implies the uniqueness. D 

As for Proposition 3.2, we only need to let Ai be those w(0) for all the solutions 
u constructed above, satisfying (3.11) with Eq = v 2 jC. 

4. One-pass theorem 

The key observation in our analysis of global dynamics is that any solution with 
energy only slightly higher than the ground states can come close to the ground 
states at most once. More precisely, if a solution u passes in and out of a small 
neighborhood of {±Q}, then it can never come back again. In particular, there is no 
homoclinic orbit connecting ±Q with themselves. Even though an orbit connecting 
Q and — Q should be called heteroclinic, we are regarding it as homoclinic, by 
identifying ±Q as one point, since there will be no difference in precluding them 
by our argument. This point of view becomes more natural if one considers NLS, 
where the ground state is really a connected set (topologically a cylinder in the 
radial case) generated from Q by the invariant group action (the modulation and 
the scaling symmetries). 

The key ideas to preclude homoclinic orbits are: 

• When a solution comes very close to Q, but does not fall on the center-stable 4 
manifold, then it eventually leaves any small neighborhood of Q in such a 
way that the unstable, i.e., exponentially growing, mode dominates all the 
others (the stable and dispersive 5 ones). This also applies to the negative 
time direction. 

• For the solutions ejected from a small neighborhood of Q with a dominating 
velocity in the unstable direction, we can use the virial identity, after suitable 
localization, as a Lyapunov-type quantity. 



If one wishes not make any reference to Proposition 3.1, then the dichotomy expressed in this 
idea simply becomes the general "trapped by Q" or "non-trapped by Q" distinction, cf. Theo- 
rem 1.1. 

However, we do not control quantitatively what "eventually" means, nor do we need to. The 
dynamics that takes place before the exponential expansion dominates is very complicated and 
relics on an interplay between the different components, and the dispersive PDE behavior can be 
of the same order of magnitude as the other dynamics. We therefore treat the "pre-exit" dynamics 
as a black box. 
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The localized virial identity will be used in the following form. For a smooth function 
w cutting-off the outside of light cones, we have 

V w (t) := (wu t \(xV + Vx)«), V w (t) = -K 2 (u(t)) + error, (4.1) 

where the error term is due to the cut-off and bounded by the linear energy at t in 
the exterior region. 

In order to use V w as a Lyapunov functional, we need a lower bound on |_ftr 2 (M)|, 
or |.Ko(w)| if we use the convexity of L 2 norm instead. In Payne-Sattinger [32] and 
Kenig- Merle [25, 26], this is achieved solely by the variational structure, which is 
not sufficient by itself in our setting. Hence, our lower bound comes in two ways: 

• While u(t) is away from but still close to ±Q, we use the hyperbolic nature 
of the eigenmode A (Lemma 4.2). 

• While u(t) is not so close or really far away from ±Q, we use the variational 
structure (Lemma 4.3). 

Note that the bound by A is not meaningful by itself once u(t) moves far away from 
±<5, since it is based on the linearization. The variational estimate is not useful 
close to ±Q, even if we knew that u(t) does not really approach ±Q, because that 
lower bound depends badly on the distance from ±Q. However, for small e > 
these two estimates exhibit sufficient overlap of their regions of validity. 

Now we introduce a nonlinear distance function to ±Q, which seems best suited 
in order to exploit the hyperbolic dynamics together with the nonlinear energy 
structure. Let 

u = a[Q + v], -y = Ap + 7, 7 _L p (4.2) 

for a = ±, and decompose the energy into the linearized part (3.7) and the higher 
order: 

E{u) - J(Q) + k 2 X 2 = \\v\\ 2 E - C(v), C(v) := (Q\v 3 ) + |M| 4 L4 /4. (4.3) 

There exists < 5e ^ 1 such that 

\\v\\ E <4S E =► |C(v)|< IHII/2. (4.4) 

Let x be a smooth function on K. such that x{ r ) = 1 f° r l r l ^ 1 an d x{ r ) = f° r 
\r\ > 2. We define 

da(u) ■= y/\\n 2 E -x(\ME/(25 E ))C(v). (4.5) 

It has the following properties 

\\v\\E/2<d a (u)<2\\v\\ E , d a (u) = \\v\\ E + 0(\\v\\ 2 E ), (4.6) 

d a (u) < S E =» d a (u) 2 = E(u) - J(Q) + k 2 X 2 . (4.7) 

Henceforth, we shall always assume that u is decomposed as in (4.2) such that 

d Q (u) := inf d±(u) = d ff (u), (4.8) 

where the choice of sign a is unique as long as dg(u) < 25 E . We also set 

\ ± (t):=\(t)±X{t)/k, (4.9) 
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the unstable/stable modes for t — > oo relative to the linearized hyperbolic evolution, 
see (3.4). First, we investigate the solutions which are close to ±Q but which are 
moving away. 

4.1. Eigenmode dominance. The first observation is that the eigenmode A be- 
comes dominant in the energy and has a fixed sign, once u is slightly away from 
the ground state, compared with the energy level. This is a static statement in the 
phase space H, and an immediate consequence of the definition of dq. 

Lemma 4.1. For any u &% satisfying 

E{u) < J{Q) + d Q (u) 2 /2, d Q {u) < S E , (4.10) 

one has dq(u) ~ |A|. 

In particular, A has a fixed sign in each connected component of the above region. 

Proof. (4.7) yields 

d Q {uf = E{u) - J{Q) + k 2 X 2 < d Q {u) 2 /2 + k 2 \ 2 . (4.11) 

and so, k 2 \ 2 /A < \\v\\ 2 E /2 < d Q {u) 2 /2 < k 2 X 2 . □ 

4.2. Ejection process. The following lemma is the key to extract the hyperbolic 
nature from our PDE. Here the linearized evolution of the eigenmode A plays the 
main role, and we specify that the solution is exiting rather than entering, by (4.13). 

Lemma 4.2. There exists a constant < Sx < o~e with the following property. Let 
u(t) be a local solution of NLKG on an interval [0,T] satisfying 

R := d Q (u(0)) < S x , E{u) < J(Q) + R 2 /2 (4.12) 

and for some t G (0, T), 

d Q (u(t))>R (0<Vt<t ). (4.13) 

Then dq(u(t)) increases monotonically until reaching Sx, and meanwhile, 
d Q {u{t)) ~ -sX(t) ~ -s\+(t) ~ e kt R, 

\X-(t)\ + \m)\\ E <R + e 2kt R 2 , (4.14) 

mm sK s (u(t)) > d Q (u(t)) - C*d Q (u(0)), 

s=0,2 

for either s = +1 or s = —1, where C* > 1 is a constant. 

Proof. Lemma 4.1 yields dq{u) — |A| as long as R < dq{u) < Se, whereas the energy 
conservation of NLKG and the equation of A give as long as dq{u) < Se, see (3.7), 

d t d Q {u) 2 = 2k 2 \\, d 2 d Q {u) 2 = 2k 2 \\\ 2 + 2k 4 \\\ 2 + 2k 2 \N p (v). (4.15) 

The exiting condition (4.13) implies d t dQ(u) 2 \ t= o > 0. Since N p (v) < \\v \\ 2 Hl , we 
have d 2 dQ(u) 2 ~ dQ(u) 2 as long as dq(u) ~ |A| <^ 1. 

Hence, imposing Sx < Se and small enough, we deduce that dq{u) > R strictly 
increases until it reaches Sx', meanwhile, dq(u) ~ sX for s G {±1} fixed. Since 

A^ - A^ = 4AA/A; > 0, (4.16) 

we also infer that A + ~ A. 
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Next, integrating the equation (3.2) for A yields 

|A(t) - A (t)| < f e k ^\N p {v{s))\ds < f e k ^\X(s)\ 2 ds, (4.17) 

Jo Jo 

where Ao denotes the linearized solution, which satisfies 

|A (t)| = |e fct A + (0) + e- fct A_(0)|/2 ~ Re kt . (4.18) 

Then by continuity in time we deduce 

X(t) ~ -sRe kt , \X(t)-X (t)\<R 2 e 2kt , (4.19) 

as long as Re kt <C 1. This yields the upper bounds on A± as well. 

To bound the remainder 7, we use the energy for the A equation, see (4.3), 

\ dt[ - k 2X/2 + \ 2 /2-C(\p)]\ = \(N p (v) -N p (\p))\\ < h\\ H i\\\ 2 . (4.20) 

Subtracting it from the energy (4.3) yields 

M*-C(t;) + C(Ap)]|<l|7ll^|A| 2 . (4.21) 

Integrating this bound and using the bound on A and 7(0), one obtains 

M-'l|2 <- d2 1 II -'II d2 2kT i A O o\ 

IM|l~£(0,T) £ R + \\l\\L?>E{0,T) R e . ( 4 - 22 ) 

which implies the desired bound on 7. 
Finally, recall from (2.14) that 

K (u) = -k 2 X(Q\p) - (2Q 3 | 7 > + 0(\\v\\ 2 Hl ), (4.23) 

and similarly we can expand K 2 around Q: 

K 2 (u) = -(k 2 /2 + 2)\{Q\p) - (2Q + g 3 | 7 ) + 0(||u||^i). (4.24) 

Since (Q\p) > by their positivity, we obtain the desired bound on K s . □ 

Note that the above proof did not use the equation for 7, see (3.2). Although 
we do have the full Strichartz estimate for the linearized evolution of 7 at our 
disposal, thanks to the spectral gap condition (1.15), the above proof does not 
require any dispersive nature of 7, and so is applicable to more general cases. This 
is indeed natural since we are dealing with that part of dynamics dominated by the 
hyperbolicity in A. 

4.3. Variational lower bounds. Now we turn to the solutions away from ±Q. The 
variational estimate is derived as an extension of Lemma 2.2. This is essentially a 
static statement, where u = (u, u) should be simply regarded as a point in the phase 
space %. 

Lemma 4.3. For any 5 > 0, there exist eo(5), k , Ki(5) > such that for any u ETi 
satisfying 

E(u) < J(Q) + e (5) 2 , d Q (u) > 5, (4.25) 

one has either 

K (u) < -Ki(6) and K 2 (u) < -Ki(6), (4.26) 
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or 

Kq(u) > mim(Ki(5) , kq\\u\\ 2 h i) and K 2 (u) > min(« 1 (5), «o||Vw|| L 2). (4.27) 

Proof, kq is an absolute constant that will be determined via the constant in (2.8). 
First we prove the conclusion separately for s = and s = 2 by contradiction. Fix 
s = or s = 2 and 5 > 0, and suppose that there exists a sequence u n G "H satisfying 
(4.25) with Eq = 1/n but neither (4.26) nor (4.27) with Ki = 1/n. Since K s (u n ) is 
bounded, the definition (2.5) implies that G s {u n ) ~ ||w n ||^-i is also bounded, and 
so K s (u n ) — > 0. Then by the same argument as in Lemma 2.1, we deduce that u n 
converges, after extraction of a subsequence, strongly to or ±Q. In the latter case, 
(4.25) implies that 

5 2 <\immf\\u n \\ 2 L 2<2el, (4.28) 

which is precluded by choosing eo(5) < 5/2. If u n — > 0, then (2.8) implies that the 
quadratic part dominates in K s (u n ) for large n, so (4.27) holds for some Kq > 
independently of 5. Thus we obtain (4.26) or (4.27), separately for s = and for 
s = 2. 

It remains to show that they have the same sign. First note that 

/C+ := {u E n | £(w) < J(Q) + 4 d Q (w) > 5, tf.(w) > 0}, 

^ (4.29) 

K~ :={ueH\ E(u) < J(Q) + eg, d (u) > 5, X s (w) < 0}, 

for s = 0, 2 are open sets satisfying 

jc+ n ^; = 0, £+ u ^o = Jc+ u ^j. (4.30) 

Since Kq and i^^ have the point in common, it suffices to show that both are 
connected. For that purpose, we use two kinds of deformations in K.+ . Fix s = 0, 2 
and take any u G /C+ satisfying 5 < oIq{u) <25<^ 5e- Recall the expansion 

2d Q {uf = k 2 X 2 + (L +7 | 7 > + Hvlli, - 2C(v), 

2{E{u) - J{Q)} = -k 2 \ 2 + (L +7 | 7 ) + ||t;||| 2 - 2C{v). ' ' ' 

Lemma 4.1 implies |A| ~ dq{u) provided that we choose eq(5) 2 < 5 2 /2. We deform 
u by increasing |A|, while fixing 7 and it. Then E{u) decreases and dq{u) increases, 
as long as dq{u) <^i Se- Meanwhile, u remains in /C+ and eventually 

E(u) - J(Q) < -d Q (u) 2 + 0(5 2 ) + o{d Q {uf) < -5 2 . (4.32) 

Thus we can deform /C+ fl {dg('u) < 26} into 

{ueH\ E(u) - J(Q) < -5 2 , K s (u) > 0} C K+, (4.33) 

which is contracted to {0} by the scaling transform as in Lemma 2.2 

u = (u, ii) ^ «, vii) {v : 1 ->■ +0). (4.34) 

For the remaining part of /C+ we also use this scaling transform, until either reaching 
0, or hitting the sphere dq(u) = 25, where it is reduced to the previous case. Thus 
we conclude that /C+ for both s are connected and coincide. □ 
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4.4. Sign function away from the ground states. The above two lemmas enable 
us to define a sign functional away from ±Q by combining those of —A and K s . 

Lemma 4.4. Let 5s '■= #x/(2C*) > where 5x and C* > 1 are constants from 
Lemma 4-2. Let < 5 < 5$ and 

H (5 ) :={uen\ E{u) < J{Q) + min(d Q (u) 2 /2, e (5) 2 )}, (4.35) 

where Sq(S) is given by Lemma 4-3- Then there exists a unique continuous function 
& : 'H(s) — > {±1} satisfying 

u G H(6), d Q {u) <5 E => 6(w) = -signA, 

u G H(s), d Q (u) > 5 => &(u) = sigaK (u) = signK 2 (u), 

where we set signO = +1 (a convention for the case u = 0). 

Proof. Lemma 4.1 implies that signA is continuous for <Iq(u) < Se, and Lemma 4.3 
implies that sign_K" (-u) = signK 2 (-u) is continuous for (Iq(u) > 8. Hence, it suffices to 
see that they coincide at (Iq{u) = 5s £ [5, Sx] in %{S)- Let u be a solution of NLKG 
with u(0) G 'His) an d d,Q(u(0)) = 5s- Then Lemma 4.2 implies that u(t) stays in 
11(8) an d signA (t) is constant, until dq{u{t)) reaches Sx, which is after — sigaK s (u(t)) 
becomes the same as signA (t), because 2C*5s < Sx- Since signi^ s (-u) is constant for 
d>Q(u) > 5, we conclude that signA (t) = — signK s (u(t)) from the beginning t — 0. D 

The & = +1 side is uniformly bounded in the energy, as the following lemma 
shows. 

Lemma 4.5. There exists M* ~ J(Q) 1//2 such that for any u E H(8 S ) satisfying 
&(u) = +1 we have \\u\\n ^ M*. 

Proof. If d,Q(u) < 5s, then 

IliTllw < HQHiri + Ml* < J{Q) 1 ' 2 + d Q {u) 1 ' 2 < JiQ) 1 ' 2 . (4.37) 

If dq{u) > 5s, then Kq(u) > together with (2.5) implies that 

\\u\\ 2 n = AE(u) - K (u) < AE{u) < 4[J(Q) + e (5 s ) 2 ] < J{Q) (4.38) 

as desired. □ 

4.5. Vanishing kinetic energy leads to scattering. Finally, we encounter the 
following problem in using the localized virial identity (4.1): in the regime K 2 (u(t)) > 
this functional can become arbitrarily small around for the u component of %. 
That is, K2(u(t)) can vanish at some time if (and only if) ||Vtt(t)||z,2 does so, see 
Lemma 4.3. Notice that we should treat this kind of vanishing only in the time aver- 
aged sense. The idea is then that all frequencies have to shift to 0, which leads to the 
scattering in both directions by the small Strichartz norm of subcritical regularity. 

Lemma 4.6. For any M > 0, there exists fio(M) > with the following property. 
Let uit) be a finite energy solution of NLKG (1.1) on [0,2] satisfying 

\\u\\ Lr (o,2;H) <M, [ ||V«(t)||| 2 dt < /i 2 (4.39) 

for some \x G (0,/xo]- Then u extends to a global solution and scatters to as 
t — > ±oo ; and moreover ||w(0I|l 3 l 6 (irxm 3 ) ^ /i 1//6 - 
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Proof. First we see that u can be approximated by the free solution 

v (t) := e i{v)t v + + e~ i{v}t v_, v± := [«(0) =F z<V) _1 u(0)]/2. (4.40) 

This follows simply from the Duhamel formula 



\V ~ u\\ L oo H i ( o,2) < ||m 3 ||l1L1(0,2) < 11^112,3^6 



(4.41) 



£ 1 1 Vm 1 1 L?Li (0,2) 1 1 Vm I k°°^(0,2) < V M < fi, 

if HqM <C 1, where we used Holder's inequality and the Sobolev embedding H 1 C L 6 . 
In particular, 

V> / ||vv(t)Hi2<ft 

Jo 

(4.42) 



= C / (CI 2 [2|t?+| 2 + 2|£_| 2 + Im{(e) _1 (e 4i ^ - 1)?+?-} 

> iivv+iiia + nvu-iii,, 



d£ 



where v denotes the Fourier transform in x of v. Now we use the Strichartz estimate 
for the free Klein-Gordon equation 

e v ' (/111 „ Ala ... *-» 



^II^bJ/^crxrs) ~ IMk> (4-43) 

where -B* denotes the Besov space with s regularity on L p . Using the terminology 

m q /L 1 I lA ZL*3 i rnpanc + LqI / _ _ 

>3' 18' 9' 



of [22, Lemma 4.1], (4.43) means that (|, ^, |) is 1-admissible. Combining it with 



Sobolev, we obtain 

i/6 (4.44) 



kllL?L|(RxR3) £ INI^/s^ 



± 

if HoM A <C 1. Therefore, we can identify u as the fixed point for the iteration in the 
global Strichartz norm 

IMlL t °°tfi(iixR3) ^ M, ||m||l3 L 6 (RxR 3) < ^ 1/6 , (4.45) 

which automatically scatters. D 

4.6. Local virial identity and non-existence of almost homoclinic orbits. 

Using the constants in Lemmas 4.2-4.6, we choose £*, 5*, R*, \x > such that 

£* <C -R* < min(5*, /t!(^) 1/2 , /% V' J(Q) 

A* < A*„(M,), / /6 « J(Q) 1/2 . (4.47) 

Suppose that a solution -u(t) on the maximal existence interval /cl satisfies for 
some £ G (0, £*], i? e (2e, i?*], and ri < r 2 < T3 E /, 

E(u) < J(Q) + e 2 , d Q (u(n)) <R< d Q {u{r 2 )) >R> d Q (u(r 3 )). (4.48) 

Then there exist T\ G (ti, 72) and T2 G (T2, T3) such that 

d Q {u(T x )) = R = d Q {u{T 2 )) < d Q {u{t)) (Ti < t < T 2 ). (4.49) 



1/2 1/2,, j/ n Ni/2x ( 4 - 46 ) 
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Lemma 4.4 gives us a fixed sign 

{±1} 3 s := e(u(t)) (T 1 <t< T 2 ). (4.50) 

Now we derive the localized virial identity with a precise error bound. The cut-off 
function is defined by 

w(tx) -h(x/(t-T 1 + S)) (t<(T 1 + T 2 )/2), 
W ^ X) -\ X (x/(T 2 -t + S)) (t>(T 1 + T 2 )/2), (4 ' 5ij 

where S ^> 1 is a constant to be determined later, and x is a radial smooth function 
on M 3 satisfying x{ x ) = 1 f° r M < 1 an d x( x ) = f° r \ x \ ^ 2. Using the equation 
we have 

V w (t) := (wu t \(xV + Vx)u), V w {t) = -K 2 (u(t)) + 0(E%(t)), (4.52) 

where Ej(t) denotes the exterior energy defined by 

E?(t) := / e°(«) tfa, e°(«) := [\ii\ 2 + \Vu\ 2 + \u\ 2 ]/2, 

JXi(t) 



x\>j[t-T 1 + S] (T 1 <t<a±Ik; 
x\>j[T 2 -t + S] (±f±<t<T 2 ), 



(4.53) 



xex.it) «=» ^ ' ' L " J ' - : ^ 2 < 



' Ti+T 2 

while the nonlinear version is denoted by 

Ej(t):= E(u)dx, E(u) := e°(u) - \u\ 4 /A. (4.54) 

We infer from the finite propagation speed that 

E 1 (t)<m a x(E 1 (T 1 ),E 1 (T 2 ))<e- 2S + ^ ||?(t)|||, (455) 

t=Ti,T 2 

where the term e -25 is dominating the tails of Q and p, due to their exponential 
decay. Hence, choosing 

5>|logi?|>l, (4.56) 

we obtain E\(t) < R 2 . Similarly, we have 

sup^(t)<i? 2 (T 1 <t<T 2 ) u 57) 

To bound the free version, we use the exterior Sobolev inequality 

IMIl4(t>s) < \\^\\m(r>s)- (4.58) 

We remark that it does not require the radial symmetry. Then 

E°(t)<Ej(t) + CE°(t) 2 (T!<t<T 2 ), (4.59) 

uniformly for j > 1. Since Ej(t) — > as j — > oo and Ej(t) < R 2 <C 1, the continuity 
in j implies that 

E°(t) ~ f?i(t) < i? 2 (T 1 <t<T 2 ). (4.60) 
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Thus we have obtained 

V w (t) = -K 2 (u(t)) + 0(R 2 ) (T!<t<T 2 ). (4.61) 

We turn to the leading term K 2 . In order to apply the ejection Lemma 4.2, 
we need the exiting property of the solution (4.13). For that purpose, take any 
t m G [Ti,T 2 ] where dq(u(t)) attains a minimum in t such that 

(R <) R m := d Q (u(t m )) = inf d Q (u(t)) < 8*. /a q 2 ) 

T\ and T 2 obviously satisfy the above, but there may be numerous other minimum 
points if u(t) is circulating in the phase space, which indeed happens for the approx- 
imating ODE obtained by the projection onto the A component. 

Applying Lemma 4.2 to u(t — t m ) and u(t m — t), as well as Lemma 4.4, one obtains 

d Q (u(t)) ~ -s\(t) ~ e^-^R^ sK 2 (u(t)) > d Q (u(t)) - C*R m , (4.63) 

until <ig («(£)) reaches Sx ^5*- Let I m denote that time interval around t m . Then 
the integral of (4.61) is estimated on each I m by 

\sV w (t)]i m > f [d Q (u(t)) ~ C*R m - 0(R 2 )]dt ~ S x , (4.64) 

thanks to the exponential growth and R < R m < 5* <C Sx- In the remainder 

I' :=[T 1 ,T 2 ]\[jl m , (4 65 ) 

m 

one has dq{u{t)) > 5*, and e < e (S*). Hence, Lemma 4.3 gives us (4.26) if s = — 1, 
or (4.27) ii s=l. 

In the latter case, Lemma 4.5 implies that \\u\\u < M* on [T±, T 2 ]. Since dg(u(t)) > 
5* ^> R for any t G I', the hyperbolic behavior (4.63) on I m implies 

[t-l,t + l]c [Ti,T 2 ], (4.66) 

and Lemma 4.6 together with (4.47) implies 

rt+l 

/ \\Vu{s)\\ 2 L 2ds> fi 2 , (4.67) 

Jt-i 

since otherwise ||it||x, 3 x, e ^ A 4 ^ 6 ^ J{Q) 1 ^ 2 , which contradicts that (Iq(u(Ti)) = 
R < J(Q) 1/2 - Combining it with (4.27) or (4.26), we obtain 

/ K 2 (u(s))ds > min(«i(^), «oA* 2 ) > #*, (4-68) 

Jt-i 

due to the choice of R in (4.46). Combining this and (4.64), we obtain 

\sV w (t)]% >S x x #{t m }. (4.69) 

The left-hand side is bounded by using the exponential decay of Q, 

< E Ht)\\Li + S\\v(t)\\ 2 E <R + SR 2 <R, (4J0) 

t=T 1 ,T 2 

provided that we choose S such that 

|logi?| <5<l/i2. (4.71) 
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Thus we have arrived at a contradiction since R <C 5x, precluding "almost homo- 
clinic" orbits, i.e., any trajectory which exits from and returns to the R neighborhood 
of {±Q}. In other words, every solution is allowed to enter and exit a sufficiently 
small neighborhood of ±Q at most once. 

Notice that the contradiction simply means that T 2 = 00, and then all the above 
analysis remains valid, except for the upper bound of [sKu]^. Thus we have proven 
the following result. 

Theorem 4.7 (One-pass theorem). Let £*, R* > be as in (4.46). If a solution u of 
NLKG on an interval I satisfies for some e G (0, £*], R G (2e, i?#], and T\ < r 2 G I, 

E{u) < J{Q) + e\ d Q (u(n)) <R = d Q (u(r 2 )), (4.72) 

then for all t G (r 2 , 00) fl / =: /', we have dq(u(t)) > R. 

Moreover, there exist disjoint subintervals ii, I 2 , • • • C I' with the following prop- 
erty: On each I m , there exists t m G I m such that 

d Q (u(t)) ~ e k ^ tml d Q (u(t m )), mmsK s (u(t)) > d Q {u{t)) - C*d Q (u(t m )), (4.73) 

s=0,2 " v ' 

where s := <B(u(t)) G {±1} is constant, dqiuit)) is increasing for t > t m , decreasing 
for t < t m , equals to 5x on dl m . For each t G I' \ \J m L m and s = 0,2, one has 
(t - 1, t + 1) C /', d Q (u(t)) > 5, and 

/•t+i 
/ miiu6r7ir a (w(0)<ft'»i^- (4-74) 

By the monotonicity, we can keep applying the above theorem at each t > t 2 until 
dg(u) reaches R*. Besides, one concludes that at any later time t m > t 2 necessarily 
rfg(-u) > R*. In other words, u cannot return to the distance R* to ±Q, after it is 
ejected to the distance 5x > R*- 

5. Blowup after ejection 

Here we prove that the solution u with &(u(t 2 )) = —1 in Theorem 4.7 blows up 
in finite time after r 2 , by the contradiction argument of Payne-Sattinger, using K . 
Suppose that u extends to all t > r 2 and let y(t) := ||-u(t)||^ 2 - From the NLKG we 
have 

y = 2[\\u\\l i+S K (u(t))]- (5.1) 

Applying the lower bound on Kq in Theorem 4.7 to the integral yields 

K>E^+ [ R 2 dt = oo, (5.2) 



1 Jr 

1 m 



and so y(t) — > 00 as t — > 00. Then from (5.1), 

y > -8E{u) + 6\\u\\l 2 + 2\\uf Hl > 6\\u\\ 2 L2 > 3(y) 2 /(2y), (5.3) 

for large t, where we used Cauchy-Schwarz for y = 2(u\u). Hence, 

d 2 t {y~ l/2 ) = -{^ /2 Y l [yy - 3(y) 2 /2] < o, (5.4) 

which contradicts that y — > oo as t — > oo. Therefore, u does not extend to t — > oo. 
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6. Scattering after ejection 

For the solution u with &{u{t2J) = +1 in Theorem 4.7, the forward global exis- 
tence follows from the energy bound Lemma 4.5. We prove its scattering to for 
t — > oo, by the contradiction argument of Kenig-Merle, using K 2 . 

Fix e G (0, £*) and let U(e,R*) be the collection of all solutions u of NLKG on 

[0, oo) satisfying 

E(u) < J(Q) + e 2 , d Q (u[0,oo)) C [i?*,oo), 6(u[0, oo)) = +1. (6.1) 

Note that the first two conditions imply that u[0, oo) C tiiSt) so that we can use 
Lemma 4.4 to define &(u). By the remark after Theorem 4.7, any solution with 
& = +1 in that theorem will eventually satisfy the above conditions. 
For each E > 0, let M(E) be a uniform Strichartz bound defined by 

M(E) := sup{|M| L 3 L 6 (0jOo) | uEU{e,R*), E(u) < E}, (6.2) 

where we chose the norm E\L\ to be an H 1 subcritical and non-sharp admissible 
Strichartz norm such that its finiteness implies scattering. We know by [22] that 
M(E) < oo for E < J(Q). In fact, in that case a uniform bound holds for Lj?L^(]R). 
In order to extend this property to J(Q) + e 2 , put 

E* = sup{E > | M(E) < oo} (6.3) 

and assume towards a contradiction that 

E*<J(Q)+e 2 . (6.4) 

We consider the nonlinear profile decomposition for any sequence u n G U(e, R*) 
satisfying 

£(«„)->£*, ||w n ||L3 L 6 (0)Oo) -)• oo. (6.5) 

We are going to show that the remainder in the decomposition is vanishing and 
there is only one profile which is a critical element, i.e., 

u±eU(e,R*), E(u i ) = E* ) ||«*||l3 L 6( 0iOO ) = oo. (6.6) 

The decomposition is given as follows. First we have the linear profile decomposition 
of Bahouri-Gerard [1]. 

Proposition 6.1. Let ip n be a sequence of free Klein-Gordon solutions bounded in 
%. Then after replacing it by a subsequence, there exist a sequence of free solutions 
yi bounded in 7i, and sequences of times iJGl such that for v 3 n and 7^ defined by 

VnW=E<W+7nW, vi(t)=V>(t + ti), (6?) 

j<k 

we have for any j < k, 7^(— t J n ) — >■ weakly in % as n — > oo ; 

lim limsup ||7nll(Lf ^nL?L«)(RxH3) = 0, lim \t 3 n - t k n \ = 00. (5.8) 
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The orthogonality of P n and the weak vanishing of 7^ implies that 



lim sup 



T, ||2 _ Vllvl^'ll 2 - ll^ll 2 
j<k 







(6.9) 



where the % norms are independent of t because all components are free solutions, 
and in particular all of v^ and 7^ are uniformly bounded in "H. It is simpler than the 
original form of the Bahouri- Gerard decomposition for the wave equation, because 
the translational symmetry does not occur here by the radial assumption, and the 
frequency parameter is fixed by the subcriticality of the remainder estimate. 

Proof. Since 7^ is bounded in H^, interpolation with the Strichartz bound implies 
that it suffices to estimate the remainder in LfL\. Let 7° := ip n and k = 0. If 

.k ._ i:„ „ ll_,fe| 



V 



limsup||7 n || L oc L 4 =0, ( 6 .io) 



n— >oo 



then we are done by putting 7^ = 7^ for all I > k. Otherwise, there exists a sequence 
tjel such that ||7*(— t*)||x,4 > u k /2 for large n. Since 7*(— t k ) G H is bounded, 
after extracting a subsequence it converges weakly in "H, and 7^(— t k ) converges 
strongly in L\. Let v k be the free solution given by the limit 

\imi k (-t k n )=v k (0), (6.11) 



.-»fe I +k\ „-?fc 



then by Sobolev \\v (0)||#i > v . We repeat the same procedure by induction 
for k = 1,2,3,.... Let U(t) denote the free Klein-Gordon propagator in %. If 
\t 3 n — t^| — )■ c G 1R for some j < k, then 

f n (-t k ) = Ufa - t k n )f n (-ti) -». 0, (6.12) 

weakly in %, hence \P n — t„\ — >■ 00 as long as v k 7^ 0. Then for all j < k, 

t +1 (-ti) = t(-ti) - v\t k n - 1{) -». (6.13) 

weakly in %. In particular we have (6.9), and so 

(6.14) 

j<k 

uniformly in k. Hence, lim sup n _ >00 ||7^||l°°l4 = v k — > 0, as k — > 00. D 



limsup \\il>n\\n > ^ \\v 3 \\ 2 H > Jj 

j<k j<k 



Before applying the profile decomposition, we translate u n in t to achieve 

d Q (u n (0)) > \&x, K 2 (u n (0)) » el (6.15) 

Since dq(u n ) remains above R*, the ejection Lemma 4.2 implies that there exists 
< T n < k~ 1 \og(5x/R*) so that d,Q(u n (T n )) > 8 X - Since & = +1, Lemma 4.5 
implies that ||w,i||l«>w(o,oo) < M*. Since ||«n||L 3 L 6 — *• °°; by the same argument as 
for (4.68), we deduce that there exists < T' n near T n such that 

K 2 {u n {T' n )) » Rl > 2el, d Q {u n {T' n )) > p x . (6.16) 

Translating u n := u n (t — T^), we obtain (6.15), in addition to (6.5). 
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Now apply the above lemma to the free solution with the same initial data as u n , 
and let w J n be the nonlinear solution with the same data as v ] n at t = 

U(t)u n (0) = J2^ + t, wi(t) = U N (t)vi(0), (617) 

j<k 

where U N (t) denotes the nonlinear Klein- Gordon propagator in %. Let P n — > P^ G 
[—00, 00] and u 3 be the nonlinear solution satisfying 

Iim||tP(0-iF'(OII« = 0, ( 6 .18) 

which exists at least locally around t — t 3 ^, as the unique solution of either the 
Cauchy problem at t 3 ^ G K. or the wave operator at t 3 ^ G {±00}. As a consequence 
of the local theory one has 

w 3 n (t) - u j (t + P n ) ^ (n ->■ 00) (6.19) 

in the energy and Strichartz norms locally around t = 0. Thus we consider the 
following nonlinear profile decomposition 

«n = 5^< + 7n + err0r > < : = «'(* + #)■ (6.20) 

Then if 2 (w n (0)) ^> e 2 , E(u n ) < J{Q) + e 2 and the orthogonality of the linear 
decomposition imply 

J{Q) -e 2 > limsup[£K) - K 2 (u n (0))/3] 

n— >oo 

>limsu P G 2 K(0)) = limsup^G 2 K(0)) + G 2 (^(0)). (6 ' 21) 

Since G 2 is positive definite, we obtain limsup n G 2 (w{(0)) < J(Q), which implies, 
via the minimizing property (2.6), that K s (u 3 n (0)) > for large n and s = 0,2. 
In particular, they all have positive energy, and so by the orthogonality applied to 
the energy, we deduce that E(u 3 ) < J(Q) except for at most one of them. Those 
nonlinear profiles scatter as t — > ±00. 

If all of the profiles scatter for t — > ±00, or more precisely if \\u 3 : \\l 3 l®(r) < °°> then 
the following long-time perturbation argument implies that the original solutions u n 
also scatter and remain bounded in the Strichartz norm uniformly in n. 

Lemma 6.2. There are continuous functions z/ , Co : (0, oo) 2 — > (0, 00) such that 
the following holds: Let I C K be an interval, u,w G C(I;H) satisfying for some 
A, B > and t G I 

\\u\\ L °°(i;H) + IH|zoo(j ;H ) < A, ||«;|| L 3 ( j. L 6) < B, (6.22) 

l|eg(w)||L t 1 (/;^) + NHIUjtfLl) + ll7o||Lf(i;L|) < MAB), (6.23) 

where eq(u) := u — Au + u — u 3 , and 70 := U(t — t )(u — w)(t Q ). Then 

||«-tw-7o||i, t ~(i;«) + ||u-H| L 3 ( j. L 6) - C (A,B)i/ . (6.24) 



(6.26) 
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Proof. Let Z := L\L X and 

l:=u -w, e : = (<9 t 2 - A + l)(u- w) -u 3 + w 3 . (6.25) 

There is a partition of the right half of /: 

t < ti < • • • < t n , ij = (tj,t j+1 ), I n (t , oo) = (t ,t n ), 

\H\ Z ( W <S (j = 0,...,n-l), n<C(B,S). 

We omit the estimate on in (— oo, t ) since it is the same by symmetry. Let 7j(t) : = 
U(t — tj) ; y(tj). Then the Strichartz estimate applied to the equations of 7 and "fj+i 
implies 

||7 - lj\\z(i 3 ) + ||7j+i - 7ilU(R) <\\(w + 7) 3 - H 3 + 4l\lI(i ,) tr n _. 

(6.27) 

< (£ + hWza,)) h\\z(i 3 ) + "o- 

Hence, by induction on j and continuity in t, one obtains 
||7lU(/ 3 ) + n + i\\z(t, +1 ,t n ) < C[\\ 7j \\ z{t . >ta) + v ] < (2C)V < (2C)V « 5, (6.28) 

provided that vq(A, B) is chosen small enough. Repeating the estimate (6.27) once 
more, we can bound the full Strichartz norms on 7. □ 

The point in using the above perturbation lemma is that the orthogonality in 
(6.8), i.e., \P n — i*| — > 00, implies that 

(E< + ^) 3 = E(<) 3 + °( 1 )' (6.29) 

j<k j<k 

in LjL 2 as n — > 00 and /c — >■ 00, as soon as each u J n has globally finite L 3 L & X norm. 
The same argument works even if the norm is finite only on (0, 00) or (— 00, 0), 
provided that the translation P n and the interval / are placed such that the infinite 
part of the norm does not contribute. Hence at least one profile does not scatter. 
Let u° be the non-scattering profile. 

If too = +00, then by the construction of the nonlinear profile, u° scatters as 
t — > 00, leading to a uniform bound on ||w n ||i, 3 L 6 (o,oo) f° r large n, a contradiction. 

If t^ = — 00, then u° scatters as t — > —00 by the same reason. Hence u° does 
not scatter as t — > +00. If <i<2(-u°(t)) > 3e on its maximal existence interval, then 
&(u°(t)) = +1 is preserved from t = +00, so u° is a global solution and a critical 
element, after time translation if necessary. Indeed, if it enters the interval (3e,R*\, 
but no deeper than that, then the ejection takes place and so the solution can never 
return to the i?* distance for later time. Suppose now that G?Q(-u°(t*)) < 3e for some 
t* G M. Then ||ii°||z, 3 z,|(— oo,t*) < °° implies that (6.29) holds in L\L? X {— 00, t* — 1°). 
Hence, the nonlinear profile decomposition becomes a good approximation of u n on 
that interval. Then taking k and n large enough, we infer from the perturbation 
lemma 

d Q {u n {t* - O) < d Q (v?(U)) + 0{e) < e < R t , (6.30) 

contradicting that d,Q(u n (0, 00)) C [i?*, 00), since t* — t^ ^ 00. 

The remaining case t^ G I is also similar. If u° scatters for t — ¥ 00, then the 
perturbation lemma implies that ||wn||L 3 L<3(o,oo) is bounded for large n. Hence, u° 
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does not scatter for t — > oo. On the other hand, the decomposition at t — gives 
\s x < \imsupd Q (u n (0)) < d Q (v°(t°J) + 0(e), (6.31) 

<J n— >oo 

and so dQ(u°(0) > 5 x /2 > 5 S > i?*. Moreover, 

K s (AO) = lim ^.(tZ°(0) > °. (6.32) 

and so we deduce from Lemma 4.4 that ©(-u ^^,)) = +1. 

By the same argument as above, m° is either a critical element or dq{ip{t*j) < 
3e for some £* > t^. Then (6.31) together with the above theorem implies that 
rfg(-u°(t)) > /?* for t < t ^. If u° does not scatter as t — > — oo, then M (t^ o — t) 
is a critical element. Suppose that u° scatters as t — > — oo. Then the profile 
decomposition is a good approximation of u n on (— oo, i* — i°) for large n. Then, by 
the perturbation lemma for large k and n, 

^ Z d Q (u°(Q) + 0(e) > limsuprf Q (M n (t* - O) > R *, (6.33) 

which contradicts e <^ R*. 

Thus in conclusion, for some T G 1R and s G {±1}, u°(st + T) is a critical element. 
Since u n is a minimizing sequence, it implies also that the other components must 
vanish strongly in the linear profile decomposition. Therefore, on a subsequence, 

lim \\u n {T' n )-v°{tl) H-H = 0, (6.34) 



n—}oo 



where T' n > is the time shift for (6.15). Both T' n and i° are bounded from above 
as n — > oo. If t° — ¥ — oo then u° scatters for t — > — oo, and so ||w n ||L 3 L 6 (-oo,T') is 
bounded for large n, by the local theory of the wave operator. 

Applying this to the sequence of solutions u n := u°(t + r n ) for arbitrary r n — >■ oo, 
one obtains the precompactness of the forward trajectory of the critical element, and 
a contradiction from a localized (time-independent) virial identity together with the 
lower bound on K 2 . These steps are simpler than those in Kenig-Merle [26] (see 
[22] for NLKG). Thus we conclude that no solution u satisfies (6.6), and therefore 
E* = J{Q)+e 2 . 

Theorem 6.3. For each e G (0, £*], there exists < M(J(Q) + e 2 ) < oo such that 
if a solution u of NLKG on [0, oo) satisfies E(u) < J(Q) + e 2 , dq(u(t)) > R* and 
&(u(t)) = +1 for all t > 0, then u scatters to as t — > oo and Hull x, 3 i 6 (o,oo) — M. 

Note that the uniform Strichartz bound is valid only on the time interval where 
the solution is already away from ±Q by a fixed distance. Without it, one cannot 
have any uniform bound even for those solutions scattering for both t — > ±oo, since 
they can stay close to ±Q for arbitrarily long time. 



7. Classification of the global behavior, proof of Theorem 1.1 

Fix < e < e* and let W := {u G H \ E(u) < J(Q) + e 2 } be the initial data set. 
We can define the following subsets according to the global behavior of the solution 
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u(t) to the NLKG equation: for a = ± respectively, 

S 6 . = {u(0) G 7i £ | u(t) scatters as at — > oo}, 

T £ = {m(0) e U £ | u(t) trapped by {±Q} for at ->■ oo}, (7.1) 

£ £ = {w(0) G % £ | «(*) blows up in at > 0}. 

The trapping for 7? can be characterized as follows, with any R G (2e, R*): 

3T > 0, Vt > T, d Q (u(t)) < R. (7.2) 

Obviously those sets are increasing in e, and have the conjugation property 

X £ T = {(«(0), -«(0)) G ?/ £ I w(0) G X^}, (7.3) 

for X = S, T, B. Moreover, S + and 7+ are forward invariant by the flow of NLKG, 
while «S_ and 71 are backward invariant. We have proven in the previous sections 
that 

H £ = S £ + UT+UB £ + =S £ _UTlUB £ _, (7.4) 

the disjoint union for each sign. It follows from the scattering theory that S± are 
open. We claim the same for B±, which is not a general fact. By the energy estimate 

\\u(t)\\ L oo H {o,T) < ||M(0)||w + r||«|||o Ofl .i (0>r) , (7.5) 

we deduce that if t — T* is the blow-up time then 

\\u(t)\\u>\T*-t\- l l\ (7.6) 

Hence, from the integral inequality 

d 2 t\\u\\% = 2[\\u\\% - K (u(t))} > 6\\ii\\ 2 Li + 2\\uf Hl - 8E(u), (7.7) 

we deduce that <9t||-u||^ 2 = 2{u\u) — > oo as well as K (u(t)) — > — oo as t — > T* — 0. 
We claim that if T** < T* is very close to blow-up time, then every solution starting 
in Bi(u(T**)) , the unit-ball around u(T**) relative to "H, necessarily blows up in the 
positive time direction. By the Payne-Sattinger argument, we only need to show 
that for any such solution, Ko(u(t)) < — K < for as long as it is defined. It is clear 
that this condition will hold initially. By Lemma 4.3 with the choice of e* in (4.46), 
we further see that it can be violated only if the solution returns to a neighborhood 
of ±Q of size 5s- However, in that case necessarily |(m|m)| < 1, which is impossible 
since (u\ii) starts off very large and has to increase as long as K (u(t)) < 0. Thus, 
any solution close to u around t = T** cannot come close to ±Q and has to blow up 
sooner or later, as claimed. Therefore, B £ ± are also open, so 7J are relatively closed 
in% £ . 

Since B £ + and S £ + are disjoint open, they are separated by 7J, that is, any two 
points from «SI and B £ + cannot be joined by a curve without passing through T+- 

In a small ball around ±Q, it is easy to see by means of the linearized flow that 
the open intersections B± fl St, B £ + fl B £ _ and S+ fl St are all non-empty for any 
e > 0. To obtain a solution u in B £ + fl St, choose initial data such that 

A(0) = 0, A(0) = k9e, 7(0) = 0, (7.8) 
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for some < $ <C 1. Then the energy constraint E(u) < J(Q) + e 2 is satisfied, and 
by the same argument as in Lemma 4.2, we have 

|A(t)-A (t)| <e 2k ^9 2 e 2 , || 7 (t)||* < Be + e^W, (7.9) 

as long as e k ^9e <C 8x, where the free solution A is given in this case 

A (£) = sinh(kt)0e. (7.10) 

Hence, u(t) exits the R* neighborhood of ±Q both in ±t > 0. Since X(t)t > for 
9s <C e k ^ e 9e <C 5x, we obtain &(u(t))t < after exiting, and so -u blows up in t > 
and scatters for t — > — oo. Therefore, £>± D St are both non-empty. 

In the same way, we construct solutions in B £ + fl B £ _ and S+ fl <S£, respectively. 
Let u^ be those solutions for which A are approximated by the free solutions 

\ ( ± ) {t) = ±cosh{kt)9e. (7.11) 

Then U W G SI n B £ _ and u<") G S% n 51. 

For e <C e fct #£ <C e 1 / 2 , the distances from the above solution u E B £ + fl <Si and 
from {±<5} are estimated by 

\\u(t) - u {+ \t)\\ E < e 2kt 9 2 e 2 + 9s « e, d Q (u(t)) > e kt 9e » e (7.12) 

Hence, there exists a curve T C "H e joining -u(t) and ir + )(i) within the region (5(w) < 
and (Iq{u) ^> e. Since -u(t) E St and t/ + )(£) G <BI, there exists a point po G 7^~nr. 
Since the solution starting from po enters the 3e-ball around ±Q as t — > — oo, and 
initially po is much further away and &(po) < 0, we conclude by the one-pass theorem 
that p E B £ + . Hence, 72 fl B £ + is non-empty as well. In the same way, we can find a 
point on the curve connecting u(t) and u^~\t) for some t < 0, which is in 7^ fl «S£. 
Therefore, 7± fl £>^ and 7± fl <Si are both not empty. Taking the limit £ — > +0, it is 
easy to observe that they contain infinitely many points on different energy levels. 
Finally, 7^ fl TI 3 Q is of course not empty, but there are infinitely many points 
besides Q, which can be seen by restricting 7i £ to the hyperplane u — 0: 

ni =0 := {(w, 0) | J(«) < J(Q) + e 2 } 3 M W(0). (7.13) 

Since -uW(O) G Bf and u( - )(0) G <S+, and there are infinitely many connecting 
curves in 7Cl =0 , there are infinitely many points in T £ D "H| =0 separating £>1 and 
•Si. The symmetry w(t) i— ¥ u(—t) implies that 

T £ nHl =0 cT £ nTl. (7.14) 

Thus we have shown that the nine solution sets are all non-empty. 

By using the scattering theory, we can prove that <Si and the energy section 

S= £ ■={ueH\ U N {t)u scatters as t -»• oo , E{u) = J(Q) + e 2 } (7.15) 

are both pathwise connected. To see this, let Uj E 5£ or S= £ for j — 0, 1 and let 
Mj + be their asymptotic profiles 

\\U N (t)uj-U(t)u j+ \\ E ^0 (t-»oo). (7.16) 

There exists a continuous curve «#+ : [0, 1] — >■ 7C such that 

\We+\\n = (1 - 0)IKo+|& + ^l|wi+llw> (7-17) 
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as well as ug solving NLKG on some (Tg, oo) such that 

\\ue(t) - U(t)ue+\\ E ->■ (t->oo). (7.18) 

Moreover, Tg and u$(Tg + l) are also continuous in 9. The scattering property implies 
that Ug(t) G S + for any t > Tg, and 

E{u e ) = \\u e+ \\ 2 H /2 = (1 - 0)£(z? o ) + OE^). (7.19) 

Let T := sup 0<6K1 Tg + 1 < oo, then -Uj and Uj{T) are connected by the flow, while 
Uo(T) and U\{T) are connected by ug(T), all included in <S+ or <S+ e . This proves the 
connectedness of those sets. 

So far, all arguments of this section have been independent of the spectral gap 
property (2.22). Now let us assume it for further investigation. Then 7J contains the 
center-stable manifold for the t > direction, as constructed in Proposition 3.1. In 
fact, it consists of the maximal backward extension of all solutions on this manifold. 
The linear approximation is the hyperplane kX = —X. Hence, the classification into 
those nine sets looks like <g> in the (A, A) phase plain around Q. Since the flow is C°° 
in both directions, 7± are indeed connected, 1-codimensional and smooth manifolds 
in the energy space. Moreover, a small ball around each point of T^ is separated by 
the hypersurface 7^ into S+ and B £ + . 

Next we prove that 7^ consists of two connected components T+<±\ defined by 

7^ (ct) = {ueW \ U N {t)u is trapped by aQ as t -»■ oo}. (7.20) 

To see that TL + \ is pathwise connected, take any couple of points and send them 
by the flow to the future where they are in 3e distance from Q. Let Uj with j — 0, 1 
be those two after translation, and let 7-,- be their asymptotic profiles 

\\U N (t)u 3 -Q-U(t)^\\ E ^0 (t^oo), E(u 3 ) = J(Q) + \\ l3 \\ 2 n /2. (7.21) 

There is a continuous 70 : [0, 1] — >■ H such that ||7e||^ = (1 — #)||7o||^ + ^ll7illw- Let 
Xg := (1 — 0)(uq — Q\p) + 0(ui — Q\p). Proposition 3.1 gives a continuous family of 
NLKG solutions u e such that (ug(0) — Q\p) = Xg and ug(t) — Q — U{t)^g — > as 
t — > 00. Then u Q and U\ are connected by the curve Ug(0) in T+, + y 

To see the separation of Tfr>) from Tf_,_^, suppose for contradiction that p is 
a common point of their boundaries. Then p G T+ £ by its relative closedness, so 
U N (t)p scatters to aQ with a — ±. Then its small neighborhood can be mapped 
by the flow into a ball around aQ. But it contains a point in Tl/ a \, contradicting 
the one-pass theorem. Thus 7^ consists of two connected components l~L±y 

We can observe also that 

dS £ + = T'U T= £ U S= £ (7.22) 

is connected. C follows from that B £ + and {E(u) > J(Q) + e 2 } are open. D is easily 
seen by perturbing the data from the right. The connectedness follows from that of 
the following three sets 

?"+(+) UT + 7 +) , S= £ UT= £ , T^UT^). (7.23) 
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Kenig-Merle's solutions [25, 26] can be classified in our notation, see (7.1), by 

S° + = S°_, B° + = B°_, T + ° = 71° = (7.24) 

(see [22] for the case of NLKG). This statement essentially follows from our theorems, 
because the energy constraint E(u) < E(Q) prohibits the sign change of K s and 
scattering to ±Q. The threshold case Theorem 1.4, i.e., the analogue of Duyckaerts- 
Merle [13, 14], requires a little more investigation in order to establish the uniqueness 
of the solutions in (5)-(9). 



Proof of Theorem 1.4- Let X 
be restated as 


± - n e >o 


xi 
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--s,b,t. 


Then Theorem 1.4 can 
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nrf = 


{±W-(s- 


-*)l 


5Gl}, 
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nr} = 


-- {±W + (s 


~t)\ 


seK}, 





for some solutions W±. To see that <S+ D B°_ = 0, let u be a solution in «S+ fl B®_ . 
Theorem 1.2 implies that there exists a finite interval I £ C H. for any e > 0, such 
that u^ 1 (B^ £ (±Q)) = I £ . Then I £ is decreasing as e — > +0, and hence by continuity 
of u, there exists iEl such that u(t) = Q, which implies u = Q, a contradiction. 

Next, any solution u in Tf corresponds to the case 7oo = in Proposition 3.1, 
and it is therefore parametrized by A(0) G KL Since X(t) — > as t — > oo, the 
uniqueness in the proposition implies that there are only three solutions converging 
to Q, modulo time translation: one with A > decreasing (for large t), another with 
A < increasing, and yet another with A = 0, i.e., Q. Let W± be the two solutions 
given by A(0) = ±i//10 and 7^ = 0. Since E(W±) = J{Q) and W± ^ Q, the sign 
of Kq is fixed on their trajectories. Indeed, Lemma 4.2 applies to them backward 
from a neighborhood of t — 00, from which we deduce that &{W±(t)) = =)=1 away 

from t — 00, and so W + G B°_ and W- 6 5". To prove that W± approach Q 
exponentially, we return to (3.5) and (3.25) with 7^ = 0, estimating 

||(A, 7 )||x* := \\e klai « {t - T > 0) \(t)\\ L i nL ~ { o,oo) + ll^^-^^WIU^oo) (7.26) 

for T — > 00. The same argument as in the proof of Proposition 3.1 yields 

ll(A,7)llx* < e- feT [|A(0)| + ||(A, 7)|||] + ||(A,7)l|x||(A, 7 )llx^ 

<e- fe ^ + HI(A,7)||x^ 

and so ||(A, 7)||x fc ~ e~ kT , as desired. D 
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Appendix A. Table of Notation 
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symbols 


description 


defined in 


HjWj'H^s) 


energy space and subsets 


(1.3), (1.9), (4.35) 


E(u),J(u) 


dynamic/static energy 


(1.2), (1.6) 


K s (u), G s {u) 


scaling derivatives of J(u) 


(2.3),(2.5) 


N(v),N p (v),N c (v) 


higher order nonlinearity around Q 


(2.12), (3.2) 


C(v) 


higher terms in energy around Q 


(4.3) 


V w (t) 


localized virial 


(4.52) 


\\v\\e 


linearized energy norm 


(3.7) 


dQ{u) 


nonlinear distance to ±Q 


(4.8) 


&{u) 


sign functional on %($) 


Lemma 4.4 


L+,uj 


linearized operator around Q 


(1.10), (3.2) 


p,k,P+ 


ground state of L + 


(1.10), (2.20) 


v,X,j,X± 


components of u around ±Q 


(4.2),(4.9) 


Se 


radius of energy nonlinearity 


(4.4) 


5x 


terminal radius of ejection 


Lemma 4.2 


O* 


constant in a bound on \K S \ 


Lemma 4.2 


e (S),Ko,Ki(S) 


constants for variational estimates 


Lemma 4.3 


Ss 


radius of variational estimates 


Lemma 4.4 


M* 


energy bound in & = +1 


Lemma 4.5 


£* 


energy threshold for one-pass theorem 


(4.46) 


<5* 


hyperbolic-variational threshold radius 


(4.46) 


R* 


exiting radius for one-pass theorem 


(4.46) 



Relation of the constants: 

2CJ S = 5 x <5 E ^:l<C*, n (M*) > (j, < J{Qf, 

£* = e (S s ) < RJ2 < mm(6 3 , «h(<$<?), «V /2 , J(Q) 1/2 ) 
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